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ABSTRACT 

Stellar evolution calculations predict the flux- weighted gravity g/T^g and ab- 
solute bolometric magnitude of blue supergiants to be strongly correlated. We use 
medium resolution multi-object spectroscopy of late B and early A supergiants 
in two spiral galaxies, NGC 300 and NGC 3621, to demonstrate the existence of 
such a relationship, which proves to be surprisingly tight. An analysis of high 
resolution spectra of blue supergiants in Local Group galaxies confirms this de- 
tection. We discuss the application of the relationship for extragalactic distance 
determinations and conservatively conclude that once properly calibrated it has 
the potential to allow for measurements of distance moduli out to 30.5 mag with 
an accuracy of 0.1 mag or better. 

Subject headings: galaxies: stellar content — galaxies: distances and redshift — 
stars: early-type 

1. Introduction 



Bright supergiants in external galaxies have been recognized as valuable extragalactic 
distance indicators for a long time, since the pioneering work of Hubble (1936). A number 
of works have attempted to calibrate photometric and/or spectroscopic signatures of blue 
supergiants (those of spectral type OBA) for this purpose, but the uncertainty in the de- 
rived distances (typically 0.4 mag or larger in distance modulus) have always been a major 
drawback for these techniques (Humphreys 1988, Tully & Wolff 1984). 

The discovery of a wind momentum-luminosity relationship (WLR) for blue supergiants 
(Kudritzki & Puis 2000) exploits the dependency of the strength of the radiation driven 
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winds of massive stars on stellar luminosity, offering a potentially more accurate distance 
indicator. While hot O stars provide so far the best calibration of the WLR (Puis et al. 
1996, Puis et al. 2002), it is the visually brightest late B and early A supergiants (My ~ —9) 
which offer the largest potential as extragalactic standard candles (Kudritzki 1998, Kudritzki 
et al. 1999). This can now be investigated for the nearby galaxies {D < 10 Mpc) with 
multiobject spectroscopy at 8-m class telescopes, as the exploratory work of Bresolin et al. 
(2001, 2002) has shown. Quantitative spectroscopy of individual BA supergiants leads to the 
determination of gravities, temperatures, metallicities and stellar wind parameters (based 
on the wind emission in Ha), which are then combined to provide distances. 

In this Letter we suggest a novel method, based on the absorption strengths of the 
higher Balmer lines formed in the photosphere. The concept is not entirely new, since a re- 
lationship between the equivalent width of H7 and My for Galactic and Magellanic Clouds 
BA supergiants, together with a temperature dependence, has already been discussed by 
several authors (Petric 1965, Crampton 1979, TuUy & Wolff 1984, Hill et al. 1986). The 
recent improvements in the modelling of A supergiant atmospheres in NLTE and the de- 
velopment of new diagnostic tools (Venn 1995a,b, 1999; Venn et al. 2001, Przybilla et al. 
2001; Przybilla & Butler 2001; Przybilla 2002) allows us to determine stellar parameters 
with unprecedented accuracy and reliability. Based on these achievements, we discuss here a 
promising application on a set of high-to- intermediate resolution spectra obtained for a sam- 
ple of A supergiants in the Milky Way, the Magellanic Clouds, NGC 300 and NGC 3621, and 
few additional objects in a handful of nearby galaxies. The theoretical concept is explained 
in §2, and we present the observational tests in §3 and §4. Future work on supergiants in 
nearby galaxies is briefly summarized in §5. 



2. Basic concept 

Massive stars during their evolution towards the Red Supergiant stage pass through 
the phase of late B and early A supergiants quickly and with roughly constant mass and 
luminosity (Meynet & Maeder 2000; Meynet et al. 1994; Heger & Langer 2000). This means 
that in this phase the stellar gravity g and effective temperature Tejf are coupled through 
the condition g/T^g = const. We call g/T^gih.e "flux-weighted gravity". Assuming that mass 
and luminosity follow the usual relation L oc M", (a ~ 3), we derive a relationship between 
absolute bolometric magnitude M(,o; and the flux- weighted gravity of the form 



-M^i^a\og{g/T^g) + h 



(1) 
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with a of the order of —3.75. This means that for these spectral types the fundamental stellar 
parameters of effective temperature and gravity are tightly coupled to the absolute magnitude 
rendering the possibility of purely spectroscopic distance determination. In the following, 
we refer to Eq. 1 as the "Flux-weighted Gravity - Luminosity Relationship (FGLR)". 

Assuming constant luminosity and, in particular, a simple (one-exponent) power law 
for the mass-luminosity relationship is, of course, a simplification. One might argue that the 
mass-loss history of supergiants and its dependence on stellar angular momentum and mctal- 
licity will complicate the situation. However, we are encouraged by detailed evolutionary 
calculations (Meynet & Macdcr 2000; Meynet ct al. 1994), which indicate for the luminosity 
and mass range of late B and early A supergiants that the amount of mass lost on the way 
from the main sequence is still relatively small and that differences in mass-loss caused by 
stellar rotation and metallicity have no substantial effects on the theoretical FGLRs derived. 

Table 1 lists the effective temperature scale for late B and early A supergiants based on 
the recent quantitative spectroscopic work cited above. We conclude that spectral classifi- 
cation allows for a temperature determination with a relative accuracy of about 4 percent. 
Then, for a given spectral type (and effective temperature) the higher Balmer lines as typical 
photospheric absorption lines allow for a rather precise determination of the gravity even for 
only medium resolution spectra of faint and distant objects (sec Fig. 1). Assuming a typical 
relative uncertainty of 0.05 for logg we estimate the uncertainty in Mhoi from Eq. (1) to be 
about 0.3 mag for a single object. This is a very encouraging estimate. Note that at this 
point we do not have to worry about absolute errors in the determination of g and Tf.ff, as 
long as we calibrate the method with a stellar sample at a known distance and apply the 
method strictly differentially. 



3. Supergiants in NGC 300 and NGC 3621 

Recently Bresolin et al. (2002) studied the population of blue supergiants in the Sculptor 
Group spiral galaxy NGC 300 at a distance of ~ 2.0 Mpc (m — M = 26.53, Frecdman ct al. 
2001). Using FORSl at the VLT medium resolution (~ 5 A) spectra of 70 blue supcrgiant 
candidates were obtained and spectral types, magnitudes and colours (the latter based on the 
work by Pietrzyhski et al. 2001) for 62 objects were presented. These observations provide 
the ideal data for a careful test of the new method. 

In our analysis we restrict ourselves to objects where we can be sure that the higher 
Balmer lines (H7, US, etc.) are not contaminated by HII region emission, i.e. we avoid 
objects with clear indication of nebular Ha emission or with a hint of nebular emission at 
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H7 above or below the 2-D stellar spectrum. In addition, we select only spectral types within 
a narrow range between B8 and A4, where we know from our recent work (Przybilla et al. 
2001; Przybilla & Butler 2001; Przybilla 2002) that our model atmosphere analysis tools 
are very reliable, in particular with regard to the relative accuracy of a strictly differential 
study. For the given spectral types wc adopt effective temperatures according to Table 1 and 
determine gravities from the higher Balmer lines as displayed in Fig. 1. We then calculate 
intrinsic colours with our model atmosphere code to determine reddening and extinction 
and use the calculated bolometric correction and the distance modulus to obtain bolometric 
magnitudes. 

The data set for NGC 300 is not the only one available to us. In a similar way, Bresolin 
et al. 2001 have used FORSl at the VLT to study 17 objects in the spiral galaxy NGC 3621 
at a distance of 6.7 Mpc {m — M — 29.08, Freedman et al. 2001). Applying the same selection 
criteria as above we can add four more objects to the sample and apply the same spectral 
analysis. 

The result of the test is displayed in Fig. 2, which shows a surprisingly tight correlation, 
as predicted by Eq.(l). The hnear regression coefficients are a — —3.85 and b — 13.73, 
the standard deviation of the the residual bolometric magnitude from this regression being 
a = 0.26 mag. The objects in NGC 3621 seem to indicate a somewhat smaller distance 
modulus (by 0.2 mag) than adopted. However, we prefer to wait for forthcoming stellar 
photometry of both NGC 300 and NGC 3621 with the Advanced Camera on board of HST, 
before we follow up on the relative distance of these two galaxies. 

4. Objects from Local Group galaxies 

The small standard deviation obtained in Fig. 2 might be an artifact resulting from the 
relatively low number of objects studied. Wc, therefore, analyze additional high resolution 
spectra presently available to us of a larger sample of objects in the Milky Way, the Magellanic 
Clouds, M31, M33 and NGC 6822. We apply the same technique as before, except for the 
three objects in the SMC, NGC 6822 and M33, which are extremely metal poor. For those 
(following Venn 1999) we do not rely on the spectral type, but determine Tgjj from the non- 
LTE ionization equilibrium of Mg I/II and N I/II. We note that the inclusion of the galactic 
objects is somewhat problematic, as their distances are more uncertain (Kudritzki et al. 
1999; Przybilla 2002). We also note that the strictly differential character of our analysis 
using the same model atmosphere and non-LTE line formation codes and the same technique 
in fitting Balmer line profiles is essential for the internal accuracy of our method. This is why 
we do not include the results of other published quantitative studies of additional objects. 
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for which we do not have spectra at our disposal. 

Fig. 3 shows the FGLR with all objects included. The regression coefficients are very 
similar (a = —3.71, b — 13.49) and the standard deviation of the residual bolometric mag- 
nitude has only slightly increased to cr = 0.28 mag. We note that the scatter increases at 
absolute magnitudes brighter than —8 mag (see discussion below). 

Fig. 3 also includes the relationship obtained from the evolutionary calculations includ- 
ing stellar rotation and mass- loss (Meynet & Maeder 2000). While the slopes practically 
agree (note that the stellar evolution relationship shows a slight curvature reflecting the 
mild change in the mass-luminosity exponent a), there is a small off-set in \og{g/T^ by 0.07 
dex. This can be the rcsTilt of a systematic effect in the determination of gravity or in the 
temperature scale, but it can also indicate a small deficiency of the evolutionary models. In 
any case, as long as an empirical calibration of the relationship with stars at known distances 
is used, the accuracy of the distance determination will not be affected. 

5. Discussion and future work 

The results presented in the previous sections are very encouraging. The flux-weighted 
gravities of late B and early A supergiants are obviously very tightly correlated with absolute 
bolometric magnitude. The application of this relationship, once properly calibrated, for ex- 
tragalactic distance determinations is straightforward. It requires multi-colour photometry 
of galaxies containing a young stellar population to identify possible blue supergiants and 
subsequent medium resolution (~ 5 A) multi-object spectroscopy (see Bresolin et al. 2001, 
2002) to determine effective temperature and gravity directly from the spectra. The spec- 
tral analysis will also yield bolometric correction (which is small for these spectral types) 
and intrinsic colour so that an accurate correction for reddening and extinction is possible. 
Application of the FGLR will then provide the absolute bolometric magnitude, which by 
comparison with the de-reddened visual magnitude will give the distance modulus. Assum- 
ing a residual scatter of o" = 0.3 mag for the FGLR (see Fig. 2 and 3) we estimate that with 
10 supergiant stars per galaxy we can achieve an accuracy of 0.1 mag in distance modu- 
lus. We are confldent that in one night of observing time we can reach down to V — 22.5 
with the existing very efficient medium resolution multi-object spectrographs attached to 
8m-class telescopes. With objects in an absolute magnitude range between —8 and —10 mag 
the FGLR method appears to be applicable out to distance moduli olm — M — 30.5 or even 
beyond. 

The restriction to medium resolution spectroscopy in the blue spectral range provides 
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significant advantages. Most importantly, contamination from sky and HII region emission 
is by far less critical than in the red, which is needed as an additional spectral range, for 
instance, for the WLR method, which requires the measurements of Ha profiles with at 
least ~ 2 A resolution. Moreover, the amount of observing time for an accurate distance 
determination is significantly reduced, if only spectra in the blue are required. 

The accurate calibration of the FGLR (and the WLR) will become the crucial clement 
of future work, before the method can be applied seriously for extragalactic distance deter- 
minations. Local Group galaxies with well determined distances provide the ideal laboratory 
for this purpose. Multi-object spectrographs with rather high spectral resolution attached 
to 8 to lOm-class telescopes such as FLAMES (VLT) or DEIMOS (Keck 2) will allow high 
quality spectra of large candidate samples in each galaxy with a rather modest amount of 
observing time. 

Such a systematic study of hundreds of blue supergiants in Local Group galaxies will 
not only provide an accurate calibration of the FGLR. It will also enable us to investigate 
important aspects of stellar evolution, which are related. The most crucial ones concern the 
role of metallicity and stellar rotation. Investigating stellar evolutionary tracks at different 
metaUicity and with different initial rotation at the main sequence (Meynet & Maeder 2000; 
Meynet et al. 1994) we find small, but noticeable effects on the theoretical FGLR. Mass-loss 
in evolutionary stages prior to the blue supergiant phase depends on metallicity (Kudritzki 
& Puis 2000) and has a small influence on the FGLR. Rotation affects the strength of mass- 
loss and the internal mixing processes and introduces a modiflcation of the mass-luminosity 
relationship in the blue supergiant stage. The effect of rotation becomes larger at higher 
luminosities. This might be the reason for the increase of the residual scatter at luminosities 
above M^oi — —8 mentioned above. 

Another very important issue is the fraction among the sample of observed blue su- 
pergiants evolving backwards to the blue after a previous phase as red supergiants. Those 
objects are expected to have lost a significant fraction of their mass as red supergiants and 
might form an additional sequence below the observed relationship. Evolutionary calcula- 
tions indicate that the relative number of those objects might depend crucially on metallicity 
and rotation. Systematic studies of Local Group galaxies at different metallicity, as proposed 
above, will allow to investigate this problem. 

In general, the observational detection of the tight relationship between fiux-weighted 
gravity and absolute bolometric luminosity is a triumph of two classical areas of astrophysics, 
stellar evolution and stellar atmospheres. It confirms the general scenario of stellar evolution 
with mass-loss and rotation away from the main sequence and the predicted mass-luminosity 
relation. It also confirms the power and accuracy of present day spectroscopic stellar diagnos- 
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tics. It is very satisfying to see the potential of these disciphnes for a significant contribution 
to quantitative extragalactic studies. 
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Table 1. Adopted temperature scale 



Spectral type 




B8 


12000 


B9 


10500 


AO 


9500 


Al 


9250 


A2 


9000 


A3 


8500 


A4 


8350 
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Fig. 1. — Fit of the Balmer lines H^,^ and Hg to Hn of the NGC 300 AOIa supergiant C6 (see 
Bresohn et al. 2002) using atmospheric models with T^ff — 9500K and log g — 1.60 (thick hne) 
and 1.65 (thin line) (gravities given in cgs units). Note that the use of information from many 
Balmer lines enhances the accuracy of the log g determination significantly. The FWHM 
of the instrumental profile is ~ 5 A, which is larger than the intrinsic width of the Balmer 
lines and the line broadening through rotation (typically < 50 km/s for A supergiants). 
Therefore, the calculations are convolved with the instrumental profile. This explains why 
gravity effects cannot be seen in the hne wings, but only in the cores. These effects in the 
hne cores reflect the changes of the integrated absorption hne strength (equivalent width) as 
a function of gravity. 
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Fig. 2. — Absolute bolometric magnitude versus logarithm of flux-weighted gravity of B8- 
to A4- supergiants in NGC 300 and NGC 3621. Note that T^g is used in units of 10^ K. 
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Fig. 3. — Same as Fig. 2 but including objects from Local Group galaxies. The relationship 
obtained from stellar evolution models at solar metallicity including the effects of rotation 
(Meynet & Maeder 2000) is also shown and labeled with the initial ZAMS-masses of the 
corresponding stellar models. 



